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Abstract 

We study the population inversion and Q-function of a two-level atom, in- 
teracting with single-mode laser light field, in a ^-analog harmonic oscillator 
trap for increasing q. For r = .003(q = e T ) the collapses and revivals of popu- 



* email : sshelly @p ositr on . r ut gers . edu 



1 



lation inversion become well denned facilitating experimental observation but 
for large r ~ 0.1 the time dependence of population inversion is completely 
wiped out. 
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I. INTRODUCTION 



Recent development of quantum groups has resulted in the construction of q- 
analogs of quite a few quantum systems including a very popular system that is the harmonic 
oscillator. A g-analog realization of harmonic oscillator has been given by Macfarlane 0] 
and Biedenharn || . In Ref . || coherent analogs of Glauber states have also been proposed. 
It is of great interest to use these quantum analog formulations to study physical systems 
and there are situations where a description in terms of q- analog structures is apparently 
more natural. The case of an atom being cooled through interaction with the radiation offers 
an interesting case study. It has been shown || recently that the atom can be cooled down 
to energies very close to it's zero point vibrational energy. At these energies the quantum 
effects due to center of mass motion of the atom are expected to be very important [0. 
Blockley et. al. || have studied the collapses and revivals of population inversion in a single 
two-level atom interacting with a classical single mode travelling light field while constrained 
to move in a one- dimensional harmonic Oscillator trap. They also discuss the possibility of 
observing the collapses and revivals experimentally. 

The interaction of the atom with the trap potential results in energy exchange 
between the internal degrees of freedom of the atom and the center of mass motion. As 
the atom slows down, the amount of energy exchanged in each step of the cooling process 
is not expected to be constant but a variable dependent on the initial energy state of the 
atomic center of mass. Presently we study the time evolution of a two-level atom in a q- 
deformed harmonic oscillator trap, in interaction with a single mode travelling light field. In 
a g-deformed harmonic oscillator trap, the energy spacing between the trap states, as seen 
by the atom, is a function of the initial vibrational state of the atom. As such, when there 
is a lot of energy associated with center of mass motion energy loss and gain occurs in large 
energy quanta. However as the atom cools down, the energy exchange takes place in smaller 
units of energy. Besides that the relative energy separation between successive states is 
determined by the value of the deformation parameter. Blockley et al || have shown that in 
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the Lamb-Dicke regime, when only the interaction between nearest neighbours is significant 
the model is similar to the Jaynes-Cummings Model(JCM) || with trap quanta playing 
a role similar to that of the light quanta in JCM. An atom with its center of mass in a 
coherent trap state initially, shows collapses and revivals of its population inversion. Our 
object is to investigate the system response for a trap potential expected to be closer to 
experimental situation in comparoson with a harmonic oscillator trap. The energy spacings 
between energy eigen states of the g-analog harmonic oscillator are not constant but are a 
function of g-deformation. As such for a given system a suitable choice of the deformation 
parameter can take us from the classical limit where the trap states are closely spaced to 
the Lamb- Dicke regime with well spaced trap states. In our earlier attempts []10| , pJ] to 
understand the physical nature of deformation in the context of pairing of nucleons, it has 
been found that the deformation amounts to simulating the nonlinearities of the problem 
or part of the interaction not included in the hamiltionian. With these results in mind 
we expect the -atom in q deformed oscillator trap model- to be a better description of the 
physical situation involved. 



II. q- ANALOG HARMONIC OSCILLATOR TRAP 

Consider a single two level atom having atomic transition frequency uo a in a 
quantized g-analog quantum harmonic oscillator trap(g-deformed harmonic oscillator trap) 
interacting with a single mode travelling light field. The creation and annihilation operators 
for the trap quanta satisfy the following quocommutation relations, 

aa + — qa + a = q~ N ; Na + — a + N = a + ; Na — aN = —a (1) 

Here N is the number operator. The operators a and a + act in a Hilbert space with basis 
vectors \n), n = 0, 1, 2, given by, 

(o + ) n 

l»> = 7Tnnrl°> (2) 

such that N \n) = n \n) . The vacuum state is a |0) = 0. We define here [x] q as 



and the g- analog factorial [n] q \ is recursively denned by 

[0] ? ! = [l]q! = 1 and [n] q \ = [n\ q [n — l] q \. It is easily verified that 

i i 
a + \n) — [n + l) q \n + 1} ; a \n) = [n] q \n — 1) (4) 

and N is not equal to a + a. Analogous to the harmonic oscillator one may define the q- 
momentum and g-position coordinate 
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The g-analog harmonic oscillator hamiltonian is given by 
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with eigenvalues 



H q ho — -hcj(aa + + a + a) (6) 



E n = ~h(j([n + l] q +[n} q ). (7) 

For r C 1 we have for the energy spacing 

E n+1 -E n ~hco + ( n + l) 2 T 2 huj. (8) 

As the atomic state changes by absorption or emission of a light quantum there is a 
change in it's vibrational state in the trap as well. We note that the trap states are not 
evenly spaced, the energy spacing being a function of deformation. Besides that as we move 
up in the number of vibrational quanta in the states the spacing between successive states 
increases. For a given region of vibrational states a suitable choice of deformation parameter 
should correspond to the Lamb Dicke regime. 



III. DYNAMICS 



The Hamiltonian for the system that consists the atom vibrating in the trap and 
interacting with a classical single- mode light field of frequency iO\. is given by || 
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H = \hu(a + a + aa + ) + \?iAa z + ]-hVt(Fa + + F*a~) (9) 

Z Zi Z 

where A = u a — ui , is the detuning parameter and Q is the Rabi frequency of the system. 
The operator F stands for exp(ikQ q ) = exp[ie(a + + a)\. The parameter e = is a function 
of the ratio of the recoil energy of the atom E r = ^^-and the characteristic trap quantum 
energy E t = hw in the limit q — > 1. Here k is the wave vector of the light field. The second 
term in the Hamiltonian refers to the energy associated with internal degrees of freedom of 
the atom, whereas the third term is the interaction of the atom with the light field. For 
nonzero values of deformation parameter q, the kinetic energy and potential energy of the 
atom in the trap are a function of deformation. In addition the interaction of the atom with 
the light field is also deformation dependent. 

The state of the system at a time t, 

= E 9m(t) \g, m> + E e m(t) \e, m) (10) 

m m 

is a solution of the time dependent Schrodinger equation 

H*(t)=ih±9(t). (11) 

The vector, \g, m), corresponds to the atom being in its ground state with its center of mass 
in m-th trap state. The label e stands for the excited state of the atom. 

The probability amplitudes g m and e m satisfy the following set of coupled equa- 
tions 

= \g m {t) H[m + 1], + [m] q ) - A) + hi ]T e n (t) (g, m\ a'F* |e, n) (12) 
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^ = h m (t) (u([m + 1], + [m\ q ) + A) + hi J2 9n(t) (e, m\ Fa + \g, n) (13) 

In order to evaluate the matrix elements (m\ F* \n) = (n\ F \ m)* , firstly we make 
use of a special case of Baker-Hausdorff Theorem to rewrite the operator F as a product of 
operators i.e we use the equality 
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exp[ie(a+ + a)] = e -(M 3 [« + >a]) e fea+ e fca 

In doing so we have substituted for the commutation relation of the operators a and 
a + the value of the commutator given in Eq.(|T]) in the limit q — > 1. This approximation 
considerably simplifies the further evaluation of matrix elements of F using the defining 
Eqs. ([!],§) for the operators a, a + and the vectors \n) . The final expression for m < n is given 

as, 



e -\'\ 2 (ie) n - m \mYV. ™ (e) 2k (-D k \n] 



( m \F\n) = - — L^lilV v ; v tl L"J<?- ( U ) 



Depending on the difference n — m the matrix element (m \F\ n) can be real, imaginary, 
positive or negative. For the special case of n = m + 1 and to first order in e the matrix 



element reduces to (m\ F \n) = ie^J [n] q . 

In the Lamb Dicke regime the model is analogous to g-analog of Jaynes-Cumming 
Model with center of mass motion quanta playing a role similar to that of the quan- 
tized radiation field. As pointed out in ref. ||, the effective Rabi frequency for trapping 
model can be easily calculated in the lowest approximation. In a similar way, consid- 
ering only those transitions involving single trap quantum exchange we can easily cal- 
culate the analog of the effective Rabi frequency in the deformed harmonic Oscillator 
trap. In case the driving field is tuned to the first vibrational sideband, A = ±u, all 
other transitions can be neglected if these are oscillating at sufficiently high frequencies. 
For uj ^> Q, the rotating wave approximation can be used and the dynamical equations 
solved to give the eigenvalues and the effective Rabi frequency. In ref. H the effective 



Rabi frequency obtained is /i(m) = y (lu + A) 2 + Vt 2 e 2 {m + 1). For the g-deformed oscilla- 
tor trap the effective frequency , to first order in e, is a function of deformation given by 



fj,(m) = J[^[cosh(2r(m + 1)) + 1] + A] 2 + 9?e 2 [m + 1],. 
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IV. INITIAL CONDITIONS 

The initial state center of mass motion of the atom is represented by g-analog 
Glauber Coherent state(GCS), while the atom is in the ground state. The q -analog of GCS 
is written as 

1, ,2 00 Q,™ 

\a) q = exp~~ 2a £__| n ) (15) 
n=oy[n],! 

with the g-exponential defined through 



ex P? = E 



00 n 



The complex parameter a determines the average number of trap quanta associated with 
the state \ot) q . 



V. POPULATION INVERSION AND QUASI-PROBABILITIES 

In figures (1) and (2), we plot the population inversion as a function of redefined 
time parameter t(= ||) for q = e T with r taking real values 0.0, 0.002, 0.003, 0.004, 0.006, 
0.01 and 0.1 respectively. Besides that we also use the parameter values e = ^ = 50 and 
A = 4 = —50. For r = 0.0 the results obtained by Blockley et. al. || are reproduced, 
when a harmonic oscillator trap is used and mean number of trap quanta is m = 16. In the 
numerical calculation the maximum value of m in eq.flTOD is restricted to m — 32. In fig.(l) 
for r = 0.002, the first revival is seen to peak earlier than the one for r = 0.0 and a well 
defined collapse appears between the first and the second revival. For a very small value of 
t = 0.003 the collapse and revival pattern is seen to become much more pronounced (well 
defined) as compared to the undeformed case. However as the deformation is increased to 
t = 0.004 and 0.006, the pattern again becomes diffuse after the first revival. The time 
interval after which the first revival occurs is seen to become shorter as r is increased as 
seen in Fig. (2). Another interesting feature of the collapses and revivals is that as the 
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deformation increases the collapse period of the system is seen to exhibit increasing level of 
coherent population trapping. For a large deformation, that in the context of this system is 
something like q > e 01 , the time dependence of popuation inversion is considerably washed 
away For still larger values of r for example r = 0.1 in fig. (2) no time dependence of the 
system is seen anymore. Similar results are obtained when imaginary values of the parameter 
r are used. It is expected because for very small deformations the values of [x] q for q = e T 
with r real are very close to the values of [x] q for r imaginary so long as the modulus of r 
is the same. 

To obtain further insight into the energy exchange between the atomic center 
of mass motion and the internal degrees of freedom of the atom, we calculate the quasi- 
probability distributions in the a r — c*i plane. The Quasi-probability or Q-function is defined 
as Q(a) = (l/iv)(a\p re d\a) q , p re d being the density matrix reduced for degrees of freedom 
of center of mass motion. Figure (3) shows in the upperv part the time evolution of Q- 
function for zero deformation at £ = 0.0, ti = 30, £ 2 = 130 and £ 3 = 160 respectively. 
Starting with the center of mass initial state as a coherent state \a ), the initially single 
peaked Q-function splits into two peaked function counterrotating in the complex plane. At 
£3 = 160 , when two peaks collide a revival of inversion oscillation occurs. A similar plot for 
t = 0.003 is shown for £ = 0.0, £1 = 30 and £2 = 129.6 respectively. An interesting feature 
of the Q-function for r = 0.003 is the presence of as many as five peaks at £ x = 30. A 
sequenece of Q-function distributions, for £ = 5, £ = 8, £ = 9, £ = 10, £ = 15 and £ = 20 shows 
the spreading out and breaking up of the quasi-probability into as many as eight distinct 
peaks(for £ = 15. The revival of Rabi- oscillations occurs when all peaks collide together 
into a single peak at £ 2 = 129.6. We also plot in Fig. (3) the Q-function for r = 0.004 at 
£ = 10 again showing as many as seven peaks. It is apparent that when trap potential is 
q-deformed the quasi-probability distribution is no longer symmetrically distributed about 
ai = 0.0 line. This asymmetry is a manifestation of unequal energy spacings between different 
trap states. 
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VI. CONCLUSIONS 



This calculation shows that in a q-deformed oscillator trap for a suitable defor- 
mation the collapse and revival patterns become well defined. As proposed by Blockley at 
el. || the theoretical predictions in this case can be verified experimentally. From our earlier 
work [TU,rT| with physical systems we may conjecture that the deformation corresponds to 
the presence of some nonlinearitis in the trapping potential. In this case small scale non- 
linearities are expected to play a beneficial role facilitating experimental observation. But 
large nonlinearities can wash out the collapse and revival pattern completely. This result 
is very important for constructing -an atom in a trap system- with a high probability for 
observing collapses and revivals of population inversion. 

In the study of interaction of a three level atom with radiation in Ref [12|] it 
is shown that by preparing the atom in a special way the time dependent collapses and 
revivals are either greatly diminished or vanish all together. The initial preparation forces 
the atomic population to remain coherently trapped in this configuration. Apparently q- 
deformed initial state is similar to a dressed atomic states. We may conclude therefore that 
the nonlinearities of the trapping potential lead to effects similar to those obtained by setting 
initial atomic state amplitudes to selected values. 

In Ref. JTB] collapses and revivals have been studied for JCM with an intensity 
dependent coupling constant by using g-analog of harmonic oscillator to represent the radi- 
ation quanta. In their calculation the collapses and revivals are seen to become more diffuse 
as the deformation is increased. Our result is different from theirs in the sense that the 
collapses and revivals become well defined. The major difference in their calculation and 
ours is that they use the zero order hamiltonian to be a harmonic oscillator Hamiltonian 
whereas we have explicitely included the zero order spectrum of g-analog oscillator. 
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FIGURES 

FIG. 1. Population inversion versus t(= g) for r = 0.0, 0.002, 0.003 and 0.004 

FIG. 2. Same as in Fig.l for r = 0.006, 0.008,0.01, 0.1. 

FIG. 3. Quasi-probability distribution in a r — oti plane at to = 0.0, t\ = 30, t<i = 130, £3 = 160. 
For r = 0.003, t 2 = 129.6 and for r = 0.004 all the peaks correspond to t = 10. 

FIG. 4. Same as in FIG. 3 for r = 0.003 at times shown in the figure. 
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